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(54) METHODE POUR CONSTrTUER UN MODELE REPRESENTATIF D'ECOULEMENTS POLYPHASIQUES DANS 
^ DES CONDUITES DE PRODUCTION PETROLIERE. 

(57j La methode proposee a pour objet de constituer un 
ffioaele representatif des ecoulements permanents et tran- 
sitoires dans des conduites d'un melange de fluides poly- 
phasiques, qui prend en compte un ensemble de variables 
definissant les proprietes des fluides et des modes d'ecou- 
lement: & phases separees, disperses, intenmrttents, et 
aussi des dimensions et pentes des conduites d'achemine- 
ment. On determine les grandeurs caractertsant recoup- 
ment en resolvant un ensemble d' equations de transport, 
une equation de conservation de masse par constituant et 
une equation de quantite de mouvement du melange, et en 
utilisant un module hydrodynamique et une caracteristique 
thermodynamique des fluides. On forme le modele en 
considerant que le melange est sensiblement a I'equilibre a 
chaque instant et que la composition du melange polypha- 
sique est variable tout le long de la condurte. 

Application a I' etude de reseau d'acheminement d'hydro- 
carbures par exemple. 
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La presente invention concerne une methode pour constituer un modele 
representatif des ecoulements stationnaires et transitoires dans des conduites d'un 
melange de fluides polyphasiques. 

5 La modelisation realisee par la methode selon l'invention, permet de tenir 

compte des phenomenes de transfert de masse entre phases et de quantites de 
mouvement entre Ies phases du melange, a partir d'un ensemble de variables 
definissant les proprietes des fluides, leurs modes d'ecoulement et aussi des 
variations de la pente des conduites par rapport a I'horizontale. Le modele ainsi 
10 constitue permet de faciliter la conception de reseaux de transfert d'effluents 
petroliers par exemple. 

La methode selon l'invention convient notamment pour modeliser le 
comportement de melanges polyphasiques d'hydrocarbures circulant dans des 
pipelines depuis des sites d'exploitation de gisements jusqu'a des sites de chargement 
15 ou de traitement par exemple. 

II est bien connu des specialistes que les modes d'ecoulement de fluides 
polyphasiques dans des tubes sont extremement varies et complexes. Les 
ecoulements diphasiques, par exemple, peuvent etre stratifies, la phase liquide 
s'ecoulant dans la partie inferieure de la conduite, ou intermittents avec une 

20 succession de bouchons liquides et gazeux, ou bien encore disperses, le liquide etant 
entraine sous forme de fines gouttelettes. Le mode d'ecoulement varie notamment 
avec Finclinaison des conduites par rappon a I'horizontale et il depend du debit de la 
phase gazeuse, de la temperature etc. Le glissement entre les phases qui varie selon 
que Ton considere les troncons ascendants ou descendants, entraine ..des variations de 

25 pression sans qu'il y ait toujours compensation. Les caracteristiques du reseau 
d'ecoulement (dimensionnement, pression, debit de gaz etc.) doivent etre 
determinees avec soin. 

Parmi les nombreuses publications portant sur le comportement des 
ecoulements notamment diphasiques dans des conduites, on peut citer par exemple : 
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- Fabre, J., et al 1983, Intermittent gas-liquid flow in horizontal or slightly inclined 
pipes, Int. Conference on the Physical Modelling of Multi-Phase Flow, Coventry, 
England, pp 233, 254; ou 


- Fabre, J., et al 1989, Two fluid/two flow pattern model for transient gas liquid 
5 flow in pipes, Int. Conference on Multi-Phase Flow, Nice, France, pp 269, 284, 
Cranfield, BHRA. 

Une methode de modelisation existe ou Ton traite les changements de phase 
par des processus iteratifs : Tetat du melange est cense connu a priori et, si eel a 
conduit, a des incoherences, apres calcul hydrodynamique, les caiculs sont repetes 
10 avec un nouvel 6tat du melange. Cette methode est lourde et peut etre la source de 
problemes de convergence. 

Une methode de modelisation appliquee a des milieux poreux est decrite par 
exemple par Eyrnard R., Gallouet T. t 1991, Traitement des changements de phases, 
dans la modelisation de gisements petroliers. Joumees numlriques de Besangon, 23 : 
15 et 24 septembre 199 1 . 

Par le brevet US 5 550 761 du demandeur, on connait une methode de 
modelisation des €coulements polyphasiques permanents ou transitoires qui prend en 
compte un ensemble de variables definissant les proprietes des fluides et des modes 
d'ecoulement, et aussi des dimensions et pentes des conduites d'acheminement. On 
20 determine les grandeurs caracterisant l'ecoulement en resolvant un ensemble 
d'equations de transport avec une equation de conservation par phase et une equation 
de quantite de mouvement du melange, et en utilisant un modele hydrodynamique et 
une caracteristique thermodynamique des fluides. 

Pour obtenir ce modele hydrodynamique, on caracterise les regimes 
25 d'ecoulement par un parametre variant entre 0 et 1 representatif de la fraction de 
l'ecoulement qui est dans un 6tat separe (les phases sont stratifiees verticalement ou 
radialement par exemple), on determine tout regime d'ecoulement au cours de la 
resolution des Equations de transport, par une comparaison de la valeur courante de 
la fraction liquide dans les bouchons et de celle des regions a mode d'ecoulement 
30 disperse, on determine aussi la vitesse des bouchons de la phase gazeuse par rapport 
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a une vitesse critique, et on impose durant la resolution des relations de fermeture, 
des contraintes de continuite aux frontieres entre les regimes, aux fractions de 
volume gazeux ainsi qu'a la vitesse de deplacement des bouchons. 

Dans la methode precedente, on a choisi une approche "par phase" (liquide- 
5 gaz) ou La conservation de la masse est traduite par une equation de conservation par 
phase et le transfert de masse entre phases est exprime par un terme de desequilibre 
proportionnel a l'ecart entre deux valeurs de fraction massique vapeur, Tune finva^ , 

correspondant a l'equilibre, qui est foumie par la thermodynamique a composition 
globale constante, l'autre etant celle calculee en tenant compte du glissement entre 
10 phases : 

ou AGKL est un facteur dependant a priori du fluide et de la configuration de 
l'ecoulement. 

On a constate a l'usage qti'il est difficile de definir une formulation de ce 
15 terme de desequilibre, qui s'applique a toutes les situations: pentes locales des 
conduites avec des points hauts et des points bas, importants transferts de masse 
entre phases. II n'existe pas de methode fiable et robuste permettant de prendre 
correctement en compte Ie terme de desequilibre entre les phases; l'approche liquide- 
gaz ou "par phase" ne donne pas de resultats satisfaisants pour traiter les cas ou se 
20 produisent d'importants transferts entre phases. 

Pour constituer un modele representatif des ecoulements permanents et 
transitoires dans des conduites d'un melange polyphasique, qui prend en compte un 
ensemble de variables definissant les proprietes des fluides. et des modes 
d'ecoulement: a phases separees, disperses, intermittents, et aussi des dimensions et 
25 pentes des conduites d'acheminement, la methode selon Tinvention comporte 
1'utilisation d'un modele hydrodynamique de type a flux de derive et d'un modele 
thermodynamique pour definir les proprietes des constituants et la resolution d'un 
ensemble d'equations de conservation de la masse par constituant, de conservation de 
la quantite de mouvement du melange et de transfert d'energie dans le melange. 
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La methode selon l'invention est caracterisee en ce que l'on forme le modele 
en considerant que le melange est sensiblement a 1'equilibre a chaque instant et que 
la composition du melange polyphasique est variable tout le long de la conduite, la 
masse de chaque constituant du melange etant definie globalement par une equation 
5 de conservation de la masse sans consideration de son etat de phase, et en ce que Ton 
utilise par exemple un schema numerique explicite en temps de facon a faciliter la 
resolution des equations du modele. 

Le traitement de l'apparition et la disparition des phases est rendu plus simple 
par cette approche compositionnelle car la masse de chaque constituant y est 
10 considered globalement sans tenir compte de son etat de phase (monophasique ou 
polyphasique). On evite ainsi les difficultes liees a 1'approche "par phase" precedente 
ou il y a une equation de conservation par phase et ou done le nombre d'equations de 
conservation de masse varie a chaque changement d'etat des constituants selon 
qu'une phase apparait ou disparait. 

15 La methode proposee permet de traiter facilement les phenomenes 

d'apparition et de disparition de phases sans rencontrer les problemes de convergence 
des solutions qui se produisent parfois avec les methodes existantes ce qui assure la 
robustesse du code. 

Suivant un mode de mise en ceuvre, la methode comporte une resolution des 
20 equations de transfert d'energie decouplee de celle concernant la conservation de la 
masse et la quantite de mouvement. 

La methode permet une resolution des equations de transfert d'energie, 
decouplee de celle concernant la conservation de la masse et la quantite de 
mouvement et elle comporte de preference I'utilisation d'un schema numerique 

25 explicite en temps, avec cette consequence tres utile que les masses- de chacun des 
constituants sont le resultat du schema numerique sans aucun appel iteratif au 
modele hydrodynamique. Connaissant les masses des constituants et la temperature, 
le modele thermodynamique int6gre determine la pression et la composition du 
melange et notamment la fraction volumique des phases. La detection de l'apparition 

30 et de la disparition des phases est ainsi rendue plus robuste. 
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Avec ce decouplage permis par la methode, on evite ainsi d'avoir a resoudre 
simultanement le modele thermodynamique et le modele hydrodynamique. Par la, on 
evite les conflits qui peuvent naitre du fait que les solutions que fournissent 
respectivement ces modeles sont a priori d'egale pertinence, et difficiles a accorder 
5 les unes avec les autres. II en resulte que la detection des phenomenes d'apparition et 
de disparition des phases est simple et robuste. 

Suivant un mode de mise en oeuvre, la methode comporte I'assimilation d'un 
melange multi-composants tel qu'un fluide petrol ier s'ecoulant dans des conduites, a 
un melange comportant un nombre plus reduit de composants et par exemple a un 
10 melange binaire equivalent (a deux constituants) ayant sensiblement la meme 
enveloppe de phase que le melange reel, de fac/on a rendre moins couteux 
l'etablissement du modele compositionnel. 

La methode comporte avantageusement I'utilisation d'un module integre de 
determination des parametres thermodynamiques (equilibre de phase et proprietes de 
15 transport) qui donne des resultats plus representatifs que ceux tires de tables 
precalculees. 

D' autres caracteristiques et avantages de la methode selon I' invention, 
apparaitront a la lecture de la description ci-apres de modes de realisation decrits a 
titre d'exemples non limitatifs, en se referant aux dessins annexes ou : 

20 - la Fig.l presente les conditions aux limites prises en compte; 

— la Fig. 2 illustre le mode de calcul stationnaire approche dans le cas ou la 
temperature est connue; 

— la Fig.3 presente l'algorithme general du calcul stationnaire dans le cas ou le calcul 
de la temperature est demande; 

25 — la Fig.4 presente l'algorithme de determination de la pression connaissant les 
masses des constituants en utilisant un flash a pression et temperatures imposees; 

— la Fig. 5 montre ralgorithme de determination de la pression connaissant les 
masses des constituants en utilisant un flash a volumes et temperatures imposes; 
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- la Fig. 6 presente l'algorithme de calcul permettant de determiner les grandeurs 
caracterisant l'ecoulement a partir des grandeurs conservatives fournies par le 
schema numerique lors du calcul trans itoire; et 

- la Fig. 7 presente l'algorithme de resolution de Tequation de transferts thermiques. 

5 I) Inconnues et equations 

La realisation du modele n composants, p phases, comporte la resolution 
d'equations de conservation de la masse pour chacun des constituants, de 
conservation de la quantite de mouvement du melange et de l'energie du melange, 
que Ton va definir ci-apres, en notant comme suit les differents parametres : 

10 LI) Inconnues: 

Xj : fraction massique du composant i dans la phase j 

c. : fraction massique totale du composant i 

Rj : fraction volumique de la phase j 

Vj : vitesse de la phase j (m/s) 

15 P : pression (Pa) 

T : temperature (K) 

Hj : enthalpie massique de la phase j 

p. : masse volumique de la phase j (kg/m 3 ) 
T w : frottement a la paroi (Pa/m) 
20 Q w : terme d'echange de chaleur a la paroi 
9 : angle de la conduite avec rhorizontale 
V M = — =H — :vitesse barycentrique du melange (m/s) 

p = Y<pj Rj masse volumique moyenne du melange. (kg/m 3 ) 

Pj R j 

X , = : fraction massique de la phase j 

P 

25 g : acceleration de la pesanteur (mis 2 ) 

S : surface de passage du fluide (wi 2 ) 

W : variables conservatives 
F : flux du schema numerique 
Q : termes sources 
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Dans I'approche compositionnelle selon l'invention, la conservation de la 
masse est verifiee pour chacun des constituants. Le transfert de masse entre phases 
n'apparalt pas explicitement dans ces equations mais il intervient par le fait que Ton 
decrit le fluide comme un melange de composition variable le long de la conduite. 
Le melange est suppose etre a l'equilibre a chaque instant. 

1.2) Definitions : 

On dSsignera dans la suite de la description de facon abregee par : 

- "flash" un sous-programme integre de calcul des proprietes thermodynamiques 
(equilibre liquide-vapeur, composition de chacune des phases) h l'aide d'une 
equation d'etat; 

- "flash (P, T)", un "flash" realise en connaissant la composition globale du 
melange, la pression et la temperature; 

- "flash (T, V)", un "flash" realise en connaissant la composition globale du 
melange, la temperature et la masse de chacun des constituants avec une 
determination de la pression au cours du calcul, de facon que les masses soient 
verifiees; 

- "flash (P, H)", un "flash" realise en connaissant la composition globale du 
melange, la pression et l'enthalpie totale, avec une determination de la temperature 
au cours du calcul, de facon que l'enthalpie soit verifiee. 

1.3) Equations 

Les equations de conservation traitees sont les suivantes : 

- une equation de conservation de la masse pour chaque constituant i : 

- une equation de conservation de quantite de mouvement du melange : 
|-(Zp y tf y V,) + ^(l P/ K,v; + P) = r„ - (Z P/ «,)*.sin0 

- une Equation d*e*nergie du melange : 
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ou Q w est Ie terme de flux d'energie a la paroL 

Pour tenir compte de l'approche compositionnelle choisie, on impose en plus 

que : 

5 £ x ' = 1 V j 

i - l 

p 

J=l 

q .= — 

Comme conditions aux limites, on peut imposer par exemple a 1'amont les 
debits massiques de chaque constituant et la temperature, et a I'aval, la pression. 

10 Comportement thermodynamique du melange 

Par application des lois de la thermodynamique, on obtient les proprietes 
physiques du fluide necessaires au code compositionnel. 

Pour une pression, une temperature et une composition globale donnSes, la 
thermodynamique permet de connaitre les fractions massiques de chacun des 

15 composants dans chacune des phases. En supplement, elle calcule la masse 
volumique de chacune des phases presentes et elle permet de deduire les taux 
volumiques de chacune des phases dans le melange global et ainsi de detecter si le 
melange est diphasique, monophasique liquide ou monophasique gaz. Des lois 
elementaires qui seront definies plus loin permettent le calcul des proprietes de 

20 transport de chacune des phases : viscosite, conductivity thermique, chaleur 
specifique, enthalpie massique, ainsi que la tension interfaciale. 


BNSDOCID: <FR 2756044A I _l_> 


9 


2756044 


La loi thermodynamique globale (calcul de I'equilibre et des proprietes du 
melange) est notee : 


^fhernw ( ^» C \ » ■ - - C *-p+ 1 ) 


proprietes de transport 
proprietes thermiques 

Loi hydrodynamique 

Dans Ie cas d'un melange, le comportement hydrodynamique est caracterise 
par Timportance du glissement entre les phases en presence : c'est-a-dire par 1'ecart 
entre les vitesses de la phase gaz et de la phase Iiquide dans le cas d'un melange 
diphasique. Le glissement depend des proprietes thermodynamiques des fluides, de 
la fraction massique de gaz, de la vitesse moyenne du melange. II est calcule par une 
fonction appelee fonction hydrodynamique qui determine egalement la configuration 
de I 'ecoulement et les termes de frottement. Cette fonction est notee : 

ou dV^Vi-V; 

Pour constituer le modele physique convenant a un 6coulement diphasique, 
on utilise par exemple les lois de fermetures suivantes bien connues des special istes : 

- pour le frottement parietal on utilise les coefficients de frottements de type 
Churchill en ecoulement turbulent et de type Poiseuille en ecoulement laminaire ; 

- pour le frottement interfacial, on utilise une loi similaire a celle de Andritsos, N. et 
Hanratty T.J.,1987, Influence of interfacial waves in stratified gas-liquid flows, 
AiChe J., Vol. 33, pp 444 - 454 ; 

- pour le diametre des bulles, on utilise une loi du type de celle proposee par Hinze, 
J.O, 1955, Fundamentals of the hydrodynamic mechanism of splitting in 
dispersion processes, AiChe J., Voll, pp 289-295 ; 
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- pour la fraction volumique de gaz dans les bouchons Iiquide, on utilise une loi 
inspiree de Andreussi P. , Bendiksen K. T 1989, An investigation of void fraction in 
liquid slugs for horizontal and inclined gas-liquid flow, Int. J. Multiphase Flow, 
Vol. 15-2, pp 937-946. 


Le terme Q^dans l'equation des transferts thermiques correspond au terme 
d'echange du a la contribution des differents modes de transfert thermiques: la 
conduction a travers la conduite (la paroi et les isolants), la convection au sein du 
fluide et Techange entre le fluide et le milieu environnant (air, sol ou mer). Pour la 
. 10 convection au sein du fluide, on tient compte de la configuration de recoupment. 

II) Proprietes physiques des fluides 

Pour caracteriser les melanges, on classe leurs grandeurs caracteristiques en 
trois groupes : 

• Les informations resultant de l'etude de 1'equilibre entre les phases: 
15 • L'etat du melange (p-phasique, triphasique, diphasique, 


5 


Modele thermique 


monophasique gaz ou monophasique Iiquide); 


• Les fractions massiques de chacun des composants dans chacune 
des phases et le taux volumique de chaque phase dans le melange; 

• Les masses volumiques de chacune des phases. 


20 


• Les proprietes de transport (utiles pour la resolution hydrodynamique): 


• Les viscosites de chacune des phases; - 

• La tension interfaciale. 


25 


• Les proprietes utiles pour la modelisation des transferts de chaleur: 

• Les enthalpies massiques de chaque phase; 

• Les conductivity thermiques dans chaque phase; 

• Les chaleurs specifiques des phases. 


Comme outils pour la modelisation thermodynamique, on peut utiliser : 
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a) des correlations: c'est-a-dire des lois simples permettant de reprfsenter 
quantitativement les phenomena physiques dont I'avantage est d'etre tres peu 
couteuses en temps de calcul, ou bien 

b) des propri6tes calculees a partir d'un programme thermodynamique 
complet avec resolution d'une equation d'etat: soit par I'interm<§diaire d'une table 
remplie par un traitement prealable, soit par I'appel a un "flash" integre a chaque fois 
que I'on a besoin de connaltre les caracteristiques du fluide. 

Presentation des correlations 

Les proprietes physiques : masses volumiques p, viscosites, et les tensions 
interfaciales a sont calculees par les formules algebriques simples. Par exemple dans 
le cas d'un melange diphasique gaz-liquide, le gaz aura un comportement proche de 
celui d'un gaz parfait et l'on utilisera les relations suivantes : 


Norm ' 


V G = Vgn, 


'orm 


P— P 

15 Qj — Q, + - LNorm 

HL HLNorm ^ 2 


V L ~ V LNorm 


a:=Z(J Nonn 


Les valeurs p r „ ,P„.T o r v \j J 1 ~ 

fournies par l'utilisateur, pour chaque simulation, de facon a reproduire au mieux le 
20 comportement du fluide modelise. 

Coefficients d'equilibre 

On ne modelise pas directement les fractions massiques de liquide et de gaz 
de chacun des constituants mais les facteurs d'equilibre K { , K i =^- pour chacun 
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des constituants, ou P°i est la pression de saturation a temperature donnee du 
composant i. 

Apparition et disparities des phases : 

Pour determiner l'etat du melange, deux notions interviennent: la pression de 
5 bulle et la pression de rosee. A composition donnee, les pressions de bulle et de 
rosee se calculent facilement. L'etat du melange depend alors de la valeur de la 
pression: 

le melange est diphasique si P ros i C <P< Pbutu 

le melange est monophasique gaz si r < 

10 le melange est monophasique liquide si P > 

Calcul des proprietes a 1'issue du flash 

Pour le calcul des proprietes physiques, on utilise par exemple le modele 
thermodynamique de Peng-Robinson avec translation de volume, bien connu des 
specialistes. Les viscosites sont calculees par la methode dite de Lohrentz Bray 
15 Clarck, les chaleurs massiques a pression constante et enthalpies, par la methode dite 
de Passut et Danner, qui est de forme polynomiale en fonction de la temperature, et 
utilise sept coefficients, la tension interfaciale, par la methode dite des parachors, ces 
parachors et 1'exposant caracteristique &ant calcules par la methode dite de Broseta, 
toutes methodes qui sont egalement bien connues des specialistes et decrites par 
20 exemple dans les publications suivantes : 

- Broseta D., et ah, 1995, Parachors in Term of Critical Temperature, Critical 
Pressure and Acentric Factor, SPE Annual Technical Conference .and Exhibition, 
Dallas, USA, 22-25 oct 1995 ; 
^ Passut C.A., et al., 1972, 1& EC, Process des. dev., 1 J, 543 (1972) ; 
25 - Peneloux A., et ah, 1982, A consistent correction for Redlich-Kwong-Soave 
volumes, Fluid Phase Equilibria , 8 (1982) pp 7J23 Elseviers Science Publishers 
(Amsterdam) ; 
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- Peng D.Y., et al. ? 1976, A new two-constant equation of state. Ind. Eng. Chem. 
Fund. 15,59-64(1976). 

Assimilation a un melange a nombre reduit de composants 

Comme on l'a vu, une caracte>istique importante de la methode selon 
I' invention est la possibility qu'elle offre d'assimiler les melanges multi-constituants 
tels que les fluides petrol iers par exemple, a un melange d'un nombre plus reduit de 
pseudo-constituants dont les proprietes sont les plus proches possibles de eel les du 
melange reel, ceci a partir de la description detaillee de la composition d'un melange 
complexe, et par exemple a un melange binaire de deux pseudo-composants. 

Utilisation des tables de proprietes 

On utilise de preference des "flash" integres pour ce qui concerne notamment 
la determination de la fraction massique vapeur, et surtout si Ton veut traiter des 
fluides a plus de deux constituants, qui donnent des resultats bien plus representatifs 
que des tables de proprietes thennodynamiques pre-remplies a l'aide d'un programme 
de calcul sur la base de la description binaire du fluide. 

m) Determination approchee de Tetat stationnaire 

Avant toute simulation, on procede a un calcul approche de Tetat initial 
stationnaire. Ce calcul permet, en partant d'une solution proche de Tetat initial reel, 
de reduire le temps de convergence pour atteindre un etat stationnaire en accord avec 
le schema numerique. 

Pour obtenir cet etat, les equations sont resolues en ne considerant pas les 
termes en derivees temporelles et en negligeant les termes d'inertie dans Tequation de 
quantite de mouvement. Les donnees utilisees sont les conditions aux limites 
indiquees sur la Fig.l, ou T est la temperature, ql,..., qn les debits massiques des 
constituants et P est la pression. 

On resout le systeme d'equations suivant : 
p 

X PjKjVjXj = pour i = l...n 
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dP 

ax 


= T w -(j.PjRj)-gsme 


x j 


proprietes de transport 
proprietes thermiques 


~ ^ ihcrmoi^ ^ > C l '•"'^n-p+l ) 


5 Les inconnues de ce probleme ainsi pose sont : P,T,c x ,...,c /l _ p+l , Rj ^dV {j ,V } . 

La demarche suivie consiste a effectuer un premier calcul de I'aval vers 
l'amont, qui permet de determiner la pression a l'amont pour un profil de temperature 
impose. Seuls les calcul s de Thydrodynamique et de la thermodynamique sont 
effectues. En partant de I'aval (Fig.2), la pression est imposee, les debits sont connus 
10 (egaux a ceux de l'amont), la temperature est imposee (cas de profil impose, ou 
estimee dans le cas contraire). 

Si le profil de temperature doit etre calcule, ralgorithme de calcul, est : 

• Calcul de I'aval vers l'amont en ay ant estime une temperature a I'aval 
en fonction des donnees sur Tenvironnement de la conduite. 

15 • Calculs de l'amont vers I'aval effectues en estimant la pression amont 

et en effectuant les calculs thermodynamique, hydrodynamique et 
thermique. La resolution se fait par une methode de calcul approche 
de type Newton a partir de la pression amont afin de retrouver la 
pression imposee a I'aval, comme schematise sur la Fig. 3. 

20 IV) Determination des etats non-stationnaires 

Les comportements instationnaires sont provoques soit par des variations des 
conditions aux limites par rapport a un £tat initial stationnaire, soit encore par la 
geomStrie irreguliere du terrain ou de Installation qui provoque des ecoulements 
instables que les specialistes designent par "terrain slugging" et "severe slugging". 
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Ddcouplage de 1'equation de transferis thermiques : 


Comme la temperature varie beaucoup moins que la pression, la composition 
et les grandeurs hydrodynamiques, il est done possible de resoudre les echanges 
thermiques de facon decouplee du calcul de la conservation de la masse et de ia 
quantite de mouvement. 

Le systeme (quantite de mouvement, masse) peut etre resolu par un schema 
numerique connu ecrit sous forme de volumes finis, expiicite en temps par exemple, 
d'ordre 1 ou 2 en espace, tel que decrit par Roe P.L., 1980, "The use of Riemann 
problem in finite difference scheme in Lectures notes of Physics 141. 

L'equation des transferts thermiques est resolue par une methode dite des 
caracteristiques. Le decouplage entre la resolution thermique et hydrodynamique, a 
pour consequence qu'a ce stade de la resolution, la temperature est connue. 

On considere les variables conservatives suivantes : 
W Mvt =^p j R j V j \ 

Les valeurs de ces grandeurs sont fournies par le schema numerique de Roe 
ci-dessus defini. 

Le systeme ^equations a resoudre est le suivant : 
W M „=lp i R j V j \ 


x j 


propriete de transport 
proprietes thermiques 


o(y u ,x J ,r, lunm „dv ij ) = o 
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Resolution pour les bords internes du maillage : 

Le schema numerique utilise permet de decoupler la resolution 
hydrodynamique de la resolution thermodynamique. Le probleme se ramene a la 
determination des grandeurs physiques a partir des grandeurs conservatives fournies 
5 par le schema numerique par la methode decrite ci-apres. 

La connaissance des masses de chacun des constituants (\^) permet de 
determiner la concentration massique de chaque composant i : 

c< =~ 

Calcul de la pression et des proprietes thermodynamiques : 

10 On peut utiliser ici une methode standard consistant a iterer sur la pression en 

f aisant des appeis a un "flash (P, T)" jusqu'a ce que la pression conduise a des masses " 
de chacun des constituants egales a celles que le schema numerique a fournies, 
suivant rorganigramme de la Fig.4. 

On peut noter qu'un gain en temps calcul important a ete obtenu dans le cas 
15 d'un appel a un "flash (T, V) M integre, en imposant des valeurs de temperatures et de 
volumes (cf. organigramme de la Fig.5); les iterations sont faites au sein de la 
thermodynamique et evitent des calculs redondants. 

Appel a la fonction hydrodynamique 

Connaissant la pression et toutes les caracteristiques du fluide, il est possible 

W 

20 d'effectuer le calcul de la vitesse barycentrique V M = . Un appel de la 

fonction hydrodynamique ®(v M , X } , T nemu> ,dV.. ) = 0 , permet de calculer les 
vitesses des phases. 

II est interessant de souligner que la resolution du modele hydrodynamique 
est identique dans la partie transitoire et dans la partie stationnaire. Le diagramme 
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general de resolution permettant d'obtenir les grandeurs physiques a partir des 
grandeurs conservatives est montre a la Fig. 6. 

Traitement des conditions aux limites 

Dans le cas general, il y a n valeurs propres positives et une valeur propre 
negative: 

A 1 < 0 < X 2 < A 3 <...< X n 

Condition a la limite aval 

Dans le cas general, il y a done une caracteristique entrante et n 
caracteristiques sortantes. La condition a la limite traduit l*information entrante. La 
pression est imposee 

N equations de compatibilite associees aux valeurs propres positives, 
traduisent rinformation sortante : 

Pour j = l,n 

n 

La loi hydrodynamique et la loi thermodynamique doivent etre 6galement 
veriflees. 

La resolution du systeme non lin£aire ainsi obtenu permet de determiner les 
masses (Wj) et la quantite de mouvement (Wj^ vt ). 

Condition a la limite amont : 

Dans le cas general, il y a n caracteristiques entrantes et 1 caracteristique 
sortante. Les conditions a la limite, sur les debits imposes, q ]9 q 2 ,—, q n de chacun 

des composants, traduisent rinformation entrante : 

Les equations a verifier sont done les suivantes : 
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tf,-(Xp;*yVO' 5=:0 pourj=l an 

^'information sortante est traduite par l'equation de compatibilite associee a 
la valeur propre negative : 

X«^+^r^^ * (3) 

5 La loi hydrodynamique et la loi thermodynamique doivent etre egalement 

verifiees. 

La resolution du systeme non lineaire ainsi obtenu permet de determiner les 
masses (Wj) et la quantity de mouvement (WMvt)- 

Remarque generate : 

10 Les differentes etapes de resolution numerique tant sur les bords que sur les 

conditions aux limites font appel a des calculs de derivees partielles. Par souci de 
robustesse, de precision et de gain en temps calcul, la plupart des deriv6es sont 
calculees de facon analytique notamment les derivees des grandeurs 
thermodynamiques. 

1 5 Transferts thermiques transitoires 

A ce stade de la resolution, la composition du m61ange et la pression sont 
connues (resolution de la conservation de la masse et de la quantite de mouvement). 

Une methode de caracteristiques (Fig.7) resout Tequation des transferts 
thermiques ; elle fournit la valeur de l'enthalpie massique du melange. 

20 L'appel a la loi thermodynamique permet de definir la temperature : 

• La methode standard consiste a faire des appels successifs a un "flash (P, T) M 
en faisant 6voluer la temperature jusqu'a ce que l'enthalpie du melange 
calculee soit identique a celle fournie par le schema numeYique. 

• Un gain de temps de calcul peut etre obtenu en ecrivant directement un 
25 "flash (P, H)". Les iterations, dans ce cas, sont faites au sein meme du 

modele thermodynamique. Le modele thermique fournit le terme d'^change 
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thermique Q int . qui est Tune des parties du terme source, Q e nth de l'equation 
de la chaleur. 

La methode selon l'invention qui permet done de model iser la variation de la 
composition d'un melange a nombre reduit de composants dans I'espace et au cours 
du temps, a ete validee experimentalement sur des cas reels. 
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REVENDICATIONS 


1) Methode pour constituer un modele representatif des ecoulements 
permanents et transitoires dans des conduites d'un melange polyphasique, comportant 
I'utilisation d'un modele hydrodynarnique de type a flux de derive et d'un modele 
thermodynamique integre pour defmir les proprietes des constituants et la resolution 
d'un ensemble d'equations de conservation de la masse, de conservation de la 
quantite de mouvement et de transfert d'energie dans le melange, caracterisee en ce 
que Ton forme le modele en considerant que le melange est sensiblement a l'equilibre 
a chaque instant et que la composition du melange polyphasique est variable tout le 
long de la conduite, la masse de chaque constituant du melange etant definie 
globalement par une equation de conservation de la masse sans consideration de son 
etat de phase. 

2) Methode selon la revendication 1, caracterisee en ce qu'elle comporte 
I'utilisation d'un schema numerique explicate en temps de facon a faciliter la 
resolution des equations du modele. 

3) Methode selon Tune des revendications 1 ou 2, caracterise en ce qu'elle 
comporte V assimilation de melanges multi-composants a des melanges constitues 
d'un nombre plus reduit de composants equivalents. 

4) Methode selon la revendication precedente, caracterisee en ce qu'elle 
comporte V assimilation de melanges multi-composants a des melanges binaires 
equivalents. 

5) Methode selon Tune des revendications precedentes, caracterisee en ce 
qu'elle comporte une resolution des equations de transfert d'energie decouplee et 
celle concernant la conservation de la masse et la quantite de mouvement 

6) Methode selon l'une des revendications pr6cedentes, caracterisee en ce 
qu'elle comporte T utilisation d'un module integre et optimise de determination des 
parametres thermodynamiques defmissant Fequilibre des phases et les proprietes de 
transport du melange. 
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